Prospect for room temperature tunneling anisotropic magnetoresistance effect: 
density of states anisotropics in CoPt systems 
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Tunneling anisotropic magnetoresistance (TAMR) effect, discovered recently in (Ga,Mn)As fer- 
romagnetic semiconductors, arises from spin-orbit coupling and reflects the dependence of the tun- 
neling density of states in a ferromagnetic layer on orientation of the magnetic moment. Based on 
ab initio relativistic calculations of the anisotropy in the density of states we predict sizable TAMR 
effects in room-temperature metallic ferromagnets. This opens prospect for new spintronic devices 
with a simpler geometry as these do not require antiferromagnetically coupled contacts on either 
side of the tunnel junction. We focus on several model systems ranging from simple hep-Co to 
more complex ferromagnetic structures with enhanced spin-orbit coupling, namely bulk and thin 
film Llo-CoPt ordered alloys and a monatomic-Co chain at a Pt surface step edge. Reliability of 
the predicted density of states anisotropics is confirmed by comparing quantitatively our ab initio 
results for the magnetocrystalline anisotropics in these systems with experimental data. 

PACS numbers: 85.75.Mm, 75.50.Cc 
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The area of condensed matter research that aims at ex- 
ploiting synergies of magnetic and semiconducting prop- 
erties in solid state systems has served as an important 
test bed for understanding basic physics and discovering 
new applications in spintronics 0, 0]. The anomalous 
Hall effect and the tunneling anisotropic magnetoresis- 
tance (TAMR) studies are two examples of the research 
that provided main motivation for the work presented 
here. The former study demonstrated 0, 0, IE IE 
that a successful theory of the anomalous Hall effect 
in (Ga,Mn)As ferromagnetic semiconductors, based on 
spin-orbit (SO) coupling effects present in the band 
structure of Bloch states, can be directly applied to 
conventional metallic ferromagnets such as Fe, and de- 
scribe quantitatively this fundamental transport coeffi- 
cient. The TAMR effect 8] is an offspring of attempts 
to develop hybrid metal/semiconductor spin- valve de- 
vices which revealed that a spin-valve-like response can 
be achieved without the seemingly fundamental switch- 
ing sequence between parallel and antiparallel magneti- 
zations in two ferromagnetic contacts with different coer- 
civities. Instead, a single ferromagnetic material with SO 
interaction is sufficient for realizing the sensing or mem- 
ory functionality through TAMR whose phenomenology 
is even richer than that of conventional giant magnetore- 
sistance or tunneling magnetoresistance effects. For ex- 
ample, both lower and higher resistance states can be 
obtained at saturation depending on the external mag- 
netic field orientation, i.e., the TAMR device can act as 
a sensor of the absolute direction of the external field . 

The TAMR was discovered in a (Ga,Mn)As/A10 x /Au 
stack 8j and confirmed by subsequent experiments in 
(Ga,Mn)As based vertical [!| and planar [l(| tunnel de- 
vices. The former experiment underlined the importance 
of high quality interfaces and barrier materials for the 
magnitude of the effect. The lithographically defined pla- 
nar nanodevice allowed to demonstrate a direct link be- 



tween the TAMR in the tunneling nanoconstriction and 
normal anisotropic magnetoresistance in the ferromag- 
netic lead. Since the latter effect also originates from SO 
coupling and has been observed in many metallic systems 
the work suggests that the TAMR has been overlooked 
in conventional room temperature ferromagnets. To ex- 
plore this possibility we follow a theoretical strategy ap- 
plied previously to (Ga,Mn)As TAMR devices which is 
based on calculating tunneling density of states (DOS) 
anisotropics in the ferromagnet with respect to the ori- 
entation of the magnetic moment and which assumes a 
proportionality between the DOS and tunneling conduc- 
tance anisotropies. 

The paper is organized as follows: We start with a sim- 
ple hep-Co crystal as a bench mark material and then 
add to the complexities of the studied structures in order 
to enhance SO coupling related effects. The next sys- 
tem we explore is bulk Llo-CoPt alloy |lllll2j. Here Co 
produces large exchange splitting resulting in the Curie 
temperature of 750 K while the heavy elements of Pt 
substantially increase the strength of SO coupling in the 
band structure of the alloy. Effects of reduced dimen- 
sionality are explored in a thin Llo-CoPt film and, reach- 
ing the ultimate nanoscale regime, in the monatomic-Co 
chain at Pt surface step edge [Tj, 0, uZ% ■ An important 
part of the presented work is a simultaneous analysis of 
magnetocrystalline anisotropies in the studied structures. 
Given the predictive nature of our theoretical conclusions 
for tunneling DOS anisotropies, calculations of a physical 
quantity originating from the same SO coupled ab initio 
band structure and directly comparable to existing exper- 
imental data is particularly valuable. Magnetocrystalline 
anisotropies, and the derived magnetization reversals in 
external magnetic fields, are also important characteris- 
tics that define functionality of a TAMR device. 

We use the relativistic version of the full-potential lin- 
earized augmented plane- wave method (FP-LAPW) [3, 
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in which SO coupling is included in a self-consistent 
second- variational procedure |l7^ . The conventional 
(von Barth-Hedin) local spin-density approximation is 
adopted, which is expected to be valid for itinerant metal- 
lic systems. The magnetic force theorem 18;) is used to 
evaluate the DOS anisotropy and the magnetocrystalline 
anisotropy energy (MAE): starting from self-consistent 
charge and spin densities calculated for the magnetic mo- 
ment Ms aligned along one of principal axes, the Ms 
is rotated and a single energy band calculation is per- 
formed for the new orientation of Ms. The DOS and 
magnetocrystalline anisotropics result from SO coupling 
induced changes in the band eigenvalues. Importantly, 
the same set of k-points has to be used for the integra- 
tion over the Brillouin zone (BZ) for accurate evaluation 
of the DOS anisotropy and of the MAE. Furthermore, 
in order to increase the accuracy in DOS evaluation, the 
smooth Fourier interpolation scheme of Pickett et al. 0] 
is used together with linear tetrahedron method |2(j . 

We start by investigating the possibility of TAMR in 
elemental transitional metal ferromagnets. Bulk hep- 
Co is a theoretically well understood uniaxial ferromag- 
net with the MAE of 60 /ieV per unit cell and the 
easy axis of magnetization along the z-[0001] crystallo- 
graphic direction [2l|. To evaluate the DOS anisotropy, 
we first performed the self-consistent FP-LAPW calcu- 
lation of the band structure of hep-Co using experimen- 
tal lattice constant values and fixing Ms along the z- 
[0001] axis. The DOS anisotropy is obtained by rotat- 
ing Ms from the z-[0001] axis to the in-plane a;- [1000] 
axis. 1200 fc-points were used for the BZ integration. 
The integral DOSs at the Fermi energy for Ms along the 
x-[1000] axis, Ni(x), and z-[0001] axis, Ni(z), are given 
in Tab. 1. Corresponding DOS anisotropy defined as, 
\N I (x) - Ni(z)\/mm(Ni(x),N I (z)) is plotted in Fig. □ 
(see black bars). The integral DOS anisotropy ~ 0.3% is 
relatively weak, similar to the weak magnetocrystalline 
anisotropy, which is a result of a small value of SO 
coupling in elemental 3d-metals. Assuming high crys- 
talline quality of an hep-Co based TAMR device with a 
large degree of in-plane momentum conservation during 
the tunneling we can improve our estimate of the effect 
by considering the tunneling DOS at the Fermi energy 
|22, N T = 1/Qbz Jbz d 3 kv z 2 (E F )S{E(k) - E F ), where 
v z = dE(k)/dk z is the group velocity component along 
the tunneling z-direction and Ep is the Fermi energy. 
The corresponding anisotropy (see blue bar in Fig. ^ is 
a factor of 4 larger than the anisotropy in Nj which is 
a trend seen previously in the (Ga,Mn)As @- The 1.3 
% anisotropy in Nt suggests a measurable, albeit weak, 
TAMR effect even with a simple hep- Co. 

A natural recipe for enhancing the TAMR effect is by 
using ordered mixed crystals of magnetic 3d and non- 
magnetic hd transition metals. In these systems the mag- 
netic atoms can produce large exchange fields polarizing 
the neighboring heavy elements with strong SO coupling. 
Pt is particularly favorable because of its large magnetic 
susceptibility. Here we consider the Llo-CoPt alloy which 



is known to have a strong MAE with the easy-axis aligned 
perpendicular to alternating Co and Pt metal layers [ll| • 
The large anisotropy results from broken cubic symmetry 
in this layered alloy and from the presence of heavy ele- 
ments of Pt, as explained above. Our previous relativistic 
FP-LAPW calculations Q2 yield the MAE of 1.03 meV, 
in a very good agreement with experimental data |ll| . 
DOS anisotropy calculations performed here follow the 
same numerical procedures as in Ref. . Starting from 
self-consistent calculations for Ms aligned along the z- 
[001] axis, we rotate Ms to the x-[110] direction and cal- 
culate the corresponding DOS. 952 fc-points in the irre- 
ducible part of BZ (3584 fc-points in full BZ) are used 
in these numerical simulations. As expected, both the 
Nj anisotropy of 1.8% and Nt anisotropy of 4.6% (see 
Tab. 1 and Fig.^) in this mixed crystal are substantially 
larger than their hep-Co counterparts. 
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Nj + Nj = Nj ot 


Pt-Surface 


0.489 + 0.485 = 


0.974 


0.385 + 0.480 = 0.865 


Co-subsurface 


0.231 + 1.637 = 


1.868 


0.187 + 1.615 = 1.802 



TABLE I: DOS in hep-Co and Ll CoPt: Total integral DOS 
Nj (1/eV) and tunneling DOS N T (eV(a.u.) 2 ) DOS at E F for 
M s || x and z axes. For Llo-CoPt film we show also layer- 
and spin- projected DOS. The fc-space convergence achieved 
in the numerical calculations, using the smooth Fourier in- 
terpolation scheme [Tg|. gives error-bars < 1% for Ni(x) in 
hep-Co, < 2.5% for N T in hep-Co, and < 0.5% for Nj(x) and 
N T in Ll CoPt. 

So far we have demonstrated prospects of room tem- 
perature TAMR and confirmed expected trends in the 
TAMR with increasing strength of SO coupling in well 
understood conventional ferromagnets. In these bulk 
systems, ab initio relativistic calculations are known to 
provide accurate description of their properties, includ- 
ing the subtle, SO coupling related magnetocrystalline 
anisotropy. The two model systems discussed in the fol- 
lowing paragraphs allow us to take a different view at 
the DOS anisotropy which might be more closely re- 
lated to the TAMR in real tunneling structures and, in 
the latter case, explores metallic TAMR in the ultimate 
nanoscale limit. Ferromagnetic electrodes in magnetic 
tunnel junctions are grown in thin films in which sur- 
faces/interfaces effects are known to play an important 
role. To analyze these effects in our context we evalu- 
ate DOS anisotropics of a free standing Llo-CoPt film 
consisting of 5 Co and 6 Pt layers with the same inter- 
layer distances as in the bulk. 400 fc-points in the two- 
dimensional BZ were used in the numerical calculations. 
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The resulting Ni anisotropy of the film is more than a 
factor of 2 larger than in the bulk Ll -CoPt. Assuming a 
vacuum tunnel barrier and realizing that tunneling trans- 
port characteristics are dominated by properties of elec- 
tronic states near the barrier, we can use the calculated 
DOS spatially resolved to individual layers to make a re- 
fined estimate of the TAMR effect in the Llo-CoPt film. 
The partial DOS corresponding to surface Pt and subsur- 
face Co monolayers are listed in Tab 1 . The anisotropy in 
the DOS of the surface Pt layer is close to 13%. We note 
that the MAE per two-dimensional unit cell is about fac- 
tor of 5 larger than the anisotropy per three-dimensional 
unit cell in the bulk Llo-CoPt. Recalling that our CoPt 
film has 5 Co and 6 Pt layers the MAE in the two struc- 
tures are very similar. Since Curie temperatures in the 
film and bulk Llo-CoPt can also be expected to be sim- 
ilar (~ 750 K) the magnetic tunnel junction based on 
thin film Llo-CoPt is a good candidate for observing the 
TAMR effect at room temperature. 
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FIG. 1: Bulk systems: integral (filled bars) 

and tunneling (half-filled bars) DOS anisotropy, 
\N i( t)(x) - N I(T) (z)\/mm(N I{T) (x),N I{T) (z)). Ll - 
CoPt thin film: integral (filled bar) and Pt-surface 
(half-filled bar) DOS anisotropy. Monoatomic Co chain: 
\Ni(x) - Ni(z)\/mm(Ni(x),Ni(z)) for A-step (filled bar) 
and \Ni(x)-Ni(y)\/min(Ni(x),Ni(y)) for B-step (half-filled 
bar). 

Recently, Gambardella et al. 01 reported ferro- 
magnetism below 15 K in monatomic-Co chains at the 
Pt(997) surface step edge. While this system is unlikely 
to lead to room temperature spintronic applications, it 
provides a unique opportunity to study both theoretically 
and experimentally magnetic and transport anisotropics 
in monatomic ferromagnetic chains. The experiments 
[l3| revealed an unexpectedly strong MAE (~ 2.0 ± 0.2 
meV/Co atom) with the easy magnetization axis directed 
along a peculiar angle of 43° towards the Pt step edge 
and normal to the Co chain. 

We analyzed the experiment by considering two possi- 
ble Pt(lll) surface step edge geometries, the (100) mi- 
crofaceted A-step and the (111) microfaceted B-step [23j. 
In the left inset of Fig. [21 we show model supercells for the 
B-step which consist of a sub-subsurface and a subsurface 
Pt layer built of 6 rows of Pt atoms, while the surface step 



is modeled by 3 rows of Pt, one Co row, and two rows of 
empty Pt sites. Vacuum was modeled by the equivalent 
of two empty Pt layers. All interatomic distances in the 
y — z-plane were relaxed using scalar-relativistic atomic 
forces |24|. This represents an important improvement 
over previous calculations [hH ITfij which assumed values 
for pure Pt. 




FIG. 2: Monatomic-Co chain: MAE as a function of the 
magnetization angle (j>. Solid lines are fits to the numeri- 
cal points [1.145 - 2.801 cos 2 (<j> - 52°)] eV (A-step) and 
[3.530 - 4.474 cos 2 (0 - 28°)] eV (B-step) . Left inset: 
schematic crystal structure, used to represent the Co chain 
at the Pt(lll) surface B-step edge. Right inset: polar angles 
<f> and 9. 

The self-consistent relativistic FP-LAPW calculations 
were performed for magnetization aligned along the z- 
axis, using 180 /c-points in a quasi-2D BZ. The MAE is 
shown in Fig. [21 As already described in Ref. 0], the 
step-edge removes any high-symmetry directions in the 
y— z-plane. The easy axis is tilted from the z-axis towards 
the Pt step edge by 52° for the A-step and by 28° for the 
B-step, in good agreement with the experimentally ob- 
served ^3 angle of 43°. The total energy difference be- 
tween states with the magnetization along the hard and 
easy axes is « 2.8 meV/Co (A-step) and 4.5 meV/Co 
(B-step); the experimental value of ~ 2 meV/Co was 
obtained from magnetic reversal measurements at tem- 
perature 45 K. Higher experimental MAE and therefore 
even better quantitative agreement between our calcula- 
tions and experiment is expected at temperatures below 
the ferromagnetic transition temperature in this system. 

The energy differences between states with Ms along 
the x-axis (along the Co-chain) or y-axis (in-plane, per- 
pendicular to the chain), and the z-axis (out-of-plane, 
perpendicular to the chain) are shown in Tab. II. For the 
A-step, the MAE is relatively weak in the y — z-plane, 
and it becomes very strong when Ms is rotated towards 
the x-axis. For the B-step, the MAE is stronger in y — z- 
plane, and becomes weaker for x — z-plane. Furthermore, 
we can evaluate the magnetic anisotropy constants by fit- 
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Ms || x axis M$ \\ y axis Ms \\ z axis 



A-step : E x>y - E z 8.54 -0.571 

Nj{Ef) 26.225 24.963 24.719 

B-step : E x , y - E z -0.190 2.639 

NijEp) 24.731 29.313 26.606 

Anisotropy constants K± K2 K3 

relaxed A-step -1.99 2.28 1.16 

relaxed B-step -0.63 -0.71 -2.00 

unrelaxed A-step (Ref. [14]) -1.70 -0.12 -1.34 

unrelaxed B-step (Ref. [15]) -0.16 -1.06 -4.81 



TABLE II: The MAE and DOS in monatomic-Co chain: To- 
tal energy differences E x — E z , E y — E z (meV), integral DOS 
(Ni(Ef), 1/eV) for M a \\ x, y and z axes, and magnetocrys- 
talline anisotropy constants (meV). 

ting our results to the total energy angular dependence 
[ljf , [if 1 cos 29 + K 2 ( 1 - cos 29) cos 20 + K 3 sin 29 sin 0] , 
where ^1,2,3 are the uniaxial anisotropy constants, and 9 
and <f> are conventional polar angles. The calculated val- 
ues of the anisotropy constants are shown in Tab. II and 
compared with previous results assuming no structural 
relaxation 0, 0] . The Co-chain anisotropy constants, 
~2 meV/Co, set a record, exceeding the bulk and surface 
anisotropics for the conventional transitional metal ma- 
terials. Similarly, the anisotropy in Ni is large. For the 
A-step, e.g., we obtained a (6.1 %) anisotropy in the in- 
tegral DOS, indicating a sizable TAMR effect. Note that 
for the one-dimensional chain we can make an alternative 
estimate of the magnitude of the TAMR based on the cal- 
culated difference between the number of energy bands at 



Ep for the two Ms directions. We obtained 25 bands for 
Ms along the z-axis and 26 bands for Ms parallel to the 
x-axis. The corresponding 4% effect is consistent with 
the estimate based on the calculated DOS anisotropies. 
For the B-step, the change in Nj for different magnetiza- 
tion orientations is particularly large for Ms along the x- 
and y-axes. The corresponding DOS anisotropy is 18.5%. 

To conclude, our theoretical results for DOS and mag- 
netocrystallinc anisotropies in CoPt structures, based 
for both quantities on the same ab initio relativistic 
band structures and in the latter one agreeing quanti- 
tatively with available experimental data, suggest siz- 
able TAMR effects in these metal ferromagnets. While 
the anisotropies in the ferromagnetic material make the 
TAMR possible, the magnitude of the effect can be very 
sensitive to parameters of the entire tunnel device, most 
notably of the geometry and crystalline quality of the 
tunnel junction, as demonstrated in (Ga,Mn)As based 
devices The integral DOS anisotropies calcu- 

lated here for the CoPt systems are larger than their 
(Ga,Mn)As counterparts. 
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